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ABSTRACT 


Coherent Anti-Stokes Raman Spectroscopy (CARS) is a new 
type of Raman Spectroscopy. The phenomenon 1s associated with 
the nonlinear conversion of two laser beams in a medium into a 
third collimated beam at the anti-Stokes frequency of the 
medium. This analysis evaluates the performance capabilities 
of CARS under various experimental conditions for the molecules 
O.,N,,H, and CO. Laser intensity fluctuations and shot noise 
are introduced as noise sources to predict signal-to-noise 
ratios (S/N). The S/N is evaluated as a function of the par- 
tial pressure of a gas to measure the performance of CARS. An 
analysis of the uncertainty in measuring the rotational tempera- 
ture 1s presented for diatomic gases. Inverse bremsstrahlung 
in plasmas is investigated as a possible process by which the 
parameters measured by CARS might be altered. An experiment 
is discussed in which an attempt was made to remove the back- 
ground signal inherent in most CARS experiments. Partial 
interferometric cancellation of two CARS signals was observed 
by the proper positioning of two sample cells in an otherwise 


Standard CARS experiment. 
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I. INTRODUCTION 


Coherent Anti-Stokes Raman Spectroscopy (CARS) is a method 
of chemical analysis yielding molecular spectroscopic informa- 
tion by utilizing nonlinear scattering phenomena. As a part 
of the growing field of nonlinear optics, CARS has been used 
Seemeasure nonlinear coefficients in solids, liquids, 
and gases [1-14]. Much of today's efforts are focused 
on the use of CARS as a nonreactive probe in environments not 
conducive to other techniques for spectral analysis [14]. 

The intent of this thesis is to assess the usefulness of 
CARS as an analytical tool using the example of diatomic gases. 
CARS is related to the more basic phenomenon of Raman Scatter- 
ing (RS). Raman scattering involves a shift in the frequency 
of a small portion of the photons incident on a Raman active 
medium. RS has been utilized extensively to obtain molecular 
spectroscopic information. CARS yields the same type informa- 
tion as RS except that the collimated CARS signal is much 
easier to detect. 

CARS was feasible for chemical analysis only after the ad- 
vent of tunable lasers of moderately high power. The initial 
work with CARS was performed by Maker and Terhune [15] in 1963, 
but only recently has the field begun to expand rapidly, as 
evidenced by the large number of published articles. To per- 
form CARS in gases requires lasers of much higher power than in 
liquids and solids. Typically, pulsed lasers on the order of 


megawatts are used. 





The use of such high power pulsed lasers has an inherent 
problem in that the power output from pulse to pulse is not 
constant. This introduces a source of noise that hinders the 
observation of weak resonances. In addition, in a mixture of 
gases, a CARS signal from other than the gas of interest is 
present as a background signal. These two factors represent 
at present a source of interference in observing weak molecular 
resonances. 

As an analytical tool CARS can yield information about 
molecular identity, concentration and rotational temperature. 
Numerical estimates for the performance capabilities of CARS 
are based on the ability to calculate the third-order nonlinear 


3 
ae 


dielectric susceptibility, . This susceptibility is com- 


posed of a resonant term, Xp and a nonresonant term, Xu’ 
Xp 1s calculable from the physical equations that define it 
but Xvr was deduced from experimental observations. 

One source of observations for Xue relied on results from 
an experiment involving Stimulated Raman Scattering (SRS). 
Appendix A contains a development of the equations related to 
SRS and shows the relationship of the SRS gain equation to 
CARS. 

A particularly difficult environment to perform convention- 
al spectral analysis with Raman scattering 1s in an ionized 
gas called a plasma. CARS is being used in plasmas and 1s 
producing better results than conventional Raman scattering 
techniques. An analysis is made to predict the extent of in- 


teraction between a plasma and the laser used for a CARS 


experiment in the plasma. In particular, the inverse 





bremsstrahlung process for the transfer of energy is ana- 
lyzed. 

An experimental attempt was made to improve the basic 
CARS experiment by the elimination of the background signal. 
A description of the experiment and the results are included 


imiethis study. 





Gree NACTORE OF THE PROBLEM 


The intent of this study was to theoretically investigate 
the analytical capabilities of CARS. The only area investi- 
gated was concerned with diatomic gases. It was desirable to 
estimate how small a concentration of gas could be detected 
and with what accuracy the rotational temperature could be 
measured. CARS has demonstrated its usefulness as a tool for 
spectral analysis. It offers orders of magnitude greater 
sensitivity than Raman scattering for low pressure gases. [16] 

Typically, CARS experiments are performed by focusing two 
pulsed laser beams of different frequencies into a sample cell. 
The beams cross at an angle that corresponds to the momentum 
conservation requirement. For gases this angle is sufficient- 
ly small that the experiment may be performed with the beams 
collinear. One of the lasers is tunable so that the difference 
in the frequency of the two may be changed. When the differ- 
ence in the frequency becomes close to an appropriate natural 
molecular frequency, the CARS signal increases markedly in 
amplitude. 

The intensity of the CARS signal is dependent on the pro- 
duct of the square of the intensity of one beam and the first 
power of the intensity of the other. Any amplitude instability 
in the beams leads to amplitude instability in the CARS signal. 
Also, the presence of an undesirable background signal masks 


very small signals of interest. 


10 





A brief sketch of several publications in the field of 
CARS has been included in this section of the study. The 
basic equations related to CARS have also been included to 
develop the framework on which the numerical calculations 


were based. 


A. BACKGROUND 

Maker and Terhune [15], in 1963, were the first to perform 
a CARS experiment. Prior to the development of good tunable 
laser sources most of the work with CARS was developed in order 
to measure nonlinear coefficients of materials. 

Yablonovitch, Bloembergen and Wynne [17] were among the 
first to use a tunable laser to generate a continuously tunable 
difference frequency. They investigated solid state materials 
such as InSb and GaAs. Levenson, Flytzanis and Bloembergen 
[1] followed closely with an experiment in diamond to measure 


y 2? 


the variations in as a function of the difference fre- 


quency. DeMartini [7] utilized CARS with H,. gas in 1972 and 


2 
Taran [4] in 1973 reported that he had used CARS as a means of 


determining the spatial distribution of H, gas ina flame. 


2 
Taran's recent work has resulted in the spectra of several 
gases, temperature measurements and spatial distributions of 


gases in flames [14]. 


B. NONLINEAR EFFECTS IN ISOTROPIC MEDIA 
Due to the nonlinear dielectric property of materials two 
laser beams may be mixed in a medium to generate a third co- 


herent beam such that 


all 





Ww, = 2w — Ws (1) 


where 03 


pump and Stokes frequency (Ww, and w 


is the anti-Stokes frequency generated from the 
3 respectively). 
The polarization vector associated with a material may be 


expressed as a power series: 


1) = 2) 22 3) 33 
P(w) = x(a) Blt x¢5) Bly #X(a) Eqay tees Fi 
where moe is the dielectric susceptibility tensor of rank 


— 
itl and E is the applied electric field. The second order 


(2) 


term yx is equal to zero for isotropic liquids and gases 


because of inversion symmetry. The lowest order term that is 


“ho 


present in a gas is therefore the third-order nonlinear 


susceptibility. There are several physical processes that in- 


y7¥37 24H, Ft Wows, 
ci = WW. > WwW 


1 D 3° etc. The last process, ZW) - Ww 


volve the third order term, e.g., 3w 
3" is the pro- 
cess referred to as CARS. 

By assuming all fields are along one axis the components 
of the tensor, yx, and the electric field vector may be 
treated as scalars. The magnitude of the electric field may 
be expressed as: 


1(k,Z = wt) 


le,e + c.c.] (3) 


Nle 


EB. (w) = 


where Es is the amplitude, kK, is the propagation constant 
(equal to win,/c), Zz is the length along the Zz axis and 


t is time in seconds. For the CARS process the term E> in 


equation (2) includes 3B, ° ES [18]. The other terms are not 


associated with CARS. 
eZ 





)) 


The resulting third order polarization, P is given 
ly 


Piedw—Kk )Z— (2. —w ) E] 
=F [3X 7 (W370 Wy “Wy E, £5 e FERC. |2, sun 


The factor of 3 has been a point of inconsistency among au- 


thors in the past. The notation for x67? (-w 1, WW) is 
consistent with that of Bloembergen [19] for the process 
W, = W + oO + OG- (>) 


Solving Maxwell's equations for plane waves with polariza- 


E2On p') the anti-Stokes gain equation is: 
gl ge * (3y §3)) ea 3)? (6) 
dz 2cn, eZ x : 
In terms of intensity, where I. = = | a pe ee ania Gace mig 


equation (6) over a length L, the result is 


eo 





+ 
| 


2 . (AKL 
2 Sin a") 
3 [3x (3) |? 1,71 _ (7) 


a ees pu NZ 
3 140° 2 AKL 7/2 


where Ak = 2k, -k5-k,- The conversion efficiency is defined 
e 


as = and 1f Ak = 0 the conversion efficiency becomes a 


2 
maximum. Because of the small dispersion of light in gas, 


CARS may be performed with w W and w. colinear. For a 


ieee 2 3 
focused beam (assumed to be a cylinder of plane waves and 


Ak = 0) the conversion efficiency becomes: 


1S 








2 Z 2. 2 
P 4 “w 2 
= = € -(2}( u 2) ae Pi : (8) 
Z L n3c 


(3) 


The third order nonlinear dielectric susceptibility, x ; 
is composed of a frequency dependent resonant part (Xp) and a 


nearly frequency independent part, the nonresonant susceptibili- 





Cy, Xyr? 
(3 
x ) = Xp + Xyr° (e2)) 
Xp can be expressed [18] as: 
a 2Nc* do _ Yr 
Y. = SQ (10) 
aes fw - (Ww, -w ic - iy(w,-w.) ] 
2 Vv ie Ne lies? 
When WW, = WO, the resonant susceptibility becomes 
2Nc* do nh 


xX = 
R Shu, * dav? 


where N is the number density of molecules, ($5) is the 
Raman cross section for scattering, A is the difference in 
Deepability Of finding a molecule in the ground and first 
vibrational states, and y is the line width, full width at 


half-maximum intensity. 


14 








III. NUMERICAL ESTIMATION OF yi?) 


(3) 


As defined in Section II.B., yx is composed of two 
parts, the resonant and nonresonant susceptibility. MThe 
resonant part is associated with the various Raman active 
Ebequencies that characterize a molecule. The equation for 
Xp has been derived [18] from considerations of a driven 


harmonic oscillator and quantum mechanics. 1s associated 


XNR 
with the electronic structure of the molecule and is nearly 
frequency independent. 


Equation (10) defines in terms of an entire vibra- 


XR 


tional band, namely the ground state. To estimate ote 


XR 
each rotational level requires an analysis of each of the 
terms in equation (10) that are functions of ee rotational 
quantum number J. As will be shown in Section IV, the signal- 
to-noise ratio for a CARS experiment is a function of the peak 
value of Xp° Being able to predict signal-to-noise ratios is 
important when deciding which analytical technique to employ 
for a certain situation. Due to the complexity of the calcu- 
lations a computer program was utilized to predict the peak 
value of Xp for several gases. 
A. RESONANT SUSCEPTIBILITY, Xp 
1. Cross Section 

The Raman cross section for scattering is typically 


quoted as the cross section for an entire vibrational band. 


For computational purposes it was necessary to apportion this 


i> 





total value among the rotational levels. The relative popula- 
tion of a rotational level was used as the apportioning factor 


Suem that 


(23) Bg (2J3+1) TOME (12) 
oY J J Q AQ mOcat 


where 


er— 5 (2J+1) exp[-BJ (J+1)/kT], 


a 
jo 


J is the rotational quantum number, is the nuclear spin 


IJ 
weighting factor and B is the rotational constant. The values 
do 
AQ 


section were the results of experiments using a ruby laser at 


used for were from Chang [20]. His measurements of cross 
5140A. The cross section has a fourth power frequency depend- 


ence Tie thus an accurate calculation of requires that 


XR 
the data be corrected for this factor. 
2. Linewidths 
On resonance the CARS signal is proportional to the 
inverse square of the linewidth. Two general extremes of 
linewidth were used, the Doppler broadened and collision 
broadened cases. The Doppler broadened lines are narrower 
than the collision (or pressure) broadened lines. 
3. Rotational Transition Frequencies 
For molecules that exhibit Raman active vibrational 
and rotational frequencies the CARS signal increases in ampli- 
tude as the difference frequency approaches one of the molecu- 
lar resonance frequencies. The rotational frequency for a 
vibration-rotation transition of a diatomic molecule is ex- 


pressed as: 


16 





Q) rene 20) XV + 1)- a Jt el) val} 


where W is the harmonicvibrational angular frequency, Xe 
is the anharmonic constant, v and J are the quantum numbers 
of vibration and rotation respectively, and Ob is the rota- 
tion-vibration interaction constant. 
4. Nuclear Spin Weighting Factor 

Homonuclear diatomic gas molecules have a statistical 
correction factor for the population of rotational levels due 
to the spin of each of the nuclei. Alternating lines in the 
rotational spectra have intensity ratios 1:(I/I+1) where I 
is the spin of each nucleus. The strong line of D.,0. and 
No are the even numbered J levels while H.'s strong lines 


are associated with odd numbered J levels. TABLE 1 summar- 


izes the degeneracy weighting factor, Ss 


TABLE 1 


NUCLEAR SPIN WEIGHTING FACTORS (G5) 


Nuclear (I) JJ BOE aa Eoe 
Gas Spin Even # J Levels Odd # J Levels 
one aes Pn : 
ek ty i) aa ee ee 
No i ik ee 


Note that CO 1s not a homonuclear molecule and thus has no 


degeneracy factor; i.e., 5 = te aL 


5. Computer Program 


To calculate as a nearly continuous function of 


XR 
the difference frequency required a computer program. A 


17? 





Fortran program was written (See page 54) to calculate the 
maximum amplitude of Xp as a function of the difference 


frequency according to the following equation: 





4g mates ot) 7k 
es) sha? 28 8S 2 we nw)? ty wn) 
2 aan) i? olen? 


Tables 2 and 3 give the results of the computations. 


B. NONRESONANT SUSCEPTIBILITY Xn 

In 1967 Rado [21] performed a series of experiments to 
measure the magnitude of the nonresonant susceptibility, Xu! 
of several gases. The impetus for the experiments was to mea- 
Seemene Peak resonant anti-Stokes power from a hydrogen cell 
and compare that to the nonresonant anti-Stokes power generated 


in a sample cell. The ratio of powers would be equal to the 


square of the ratio of the susceptibilities: 


p x (Sample) )* 
sample , ( NR : 
Pi ( a 
2 - 


The value Rado chose to use for Xp was from an earlier 


GS) 


work with phonon lifetimes and SRS gain coefficients. The 
relationship of Xp and the SRS gain coefficient is discussed 


in Appendix A. His calculated value of Xue ECR aly is 


=o cm? ° 
Ney Sop ey) —— . Based on Table 3, x for N at 300 K 
erg 3 P 7 
is pees aaa - Thus, the ratio of the peak to nonresonant 


susceptibilities is 


18 





TABLE 2 


3 
Gas Linewidth (cm *) Temp (°K) Xp (em’/exg) 
0034 300 gcc NOR: 
0 
-16 
_ 0076 ___ 1500 ZO 
___..0055 300 gee ones 
N nia 
= 
012 1500 11 x10 
036 300 TOON 
ero 1500 4.3x10 
005 300 3.3x10 > 
CO 
07 1500 ee ie 


io 





Gas 


CO 


TABLE 3 


Linewidth (ene) 


Pa 


2007 


-14 


06 


203 


208 


ES 


Oe 


20 


1500 


1500 


(Is) 


Temp 


300 


1500 


Biwi 


1500 


300 


300 


Xp (cm? /erg) 


101° 


1971? 


1971 


1072? 


190714 


19716 


101° 


1071° 





An independent method to estimate the magnitude of Xu 
was provided by the spectra of No published by Taran [14]. 
For the unresolved spectra the difference in frequency of the 
maximum and minimum values of the CARS signal, 6éw, is given 


by the relation 





i 
Xp 2 
dw = y{l + 5 : (5G;) 
4XnR 
Assuming a linewidth of 1.8 emt and a 6w of 42 em + 
the calculated result is that 
x 
= ~ 47. (17) 
NR N 
2 


By estimating the level of the background Signal in the 


No spectrum, another ratio of Xp/Xyp 1s equal to 


Xp 
XuR 


ie 


1G (18) 


Knowing the peak value of Xp yields an estimate of XvRr- 
An earlier work by Taran [4] provides another means of 
estimating the ratio of the susceptibilities. The total 
pressure of a H,N, mixture waS maintained at one atmosphere. 
The anti-Stokes intensity generated was plotted as a function 


of the H. concentration. 


Za 


ie Xp and XyR scale linearly with pressure, then the 


total susceptibility of the mixture may be approximated by 


3 
aed 


XpP + [1-P] Xyp ee) 


where P is the partial pressure of H.- As a ratio of Xp 


to Xyp the observed anti-Stokes power is approximated by 
x 2 
(eS e2 iese Re Z 
I, a |x |" = ap | 1] P +1 ‘ (20) 


When the log of intensity is plottedagainst the log of P 
(See Figure 1), the predicted and observed curves are in fair 


agreement for a ratio of RO RreaTe: That is 


y (Hg) 
Ro ~ 10°. (21) 
(N.) 
xX 
NR 
Based on the analysis of the H,/N, plot and using the value 
of Xp om H. PeOMetaole 5, the Value Of XyR 1(G)IG N, 1s 


3 
Meeinaced to be 1.2x10~° cn, 


Thus, another estimate of the 


Gatao is 


Xp _ 2.4x1077° 


ine 
ANR Ie 250110 : 


= 200. 


The factor of four between the high and low estimates of the 
ratio is not surprising based on the variety of sources used 


to estimate them. 


Ze 








Figure 1. Taran's results follow dotted line. The solid 


‘ line is the predicted curve for -——= 10°. The 


two curves are equal from Ce to lO” ppm. 


Relative Anti-Stokes Intensity 
Arbitrary Units 


102 104 10° 


Ha Concentration (ppm) 


Be 





IV. STATISTICAL CONSIDERATIONS 


The anti-Stokes power resulting from a tightly focused 


CARS experiment is expressed as [18]: 





- (3) 2 2 
P,=K | 3x | Py P. (22) 
where 
fo 2 4m w, 2 
Ke= (2.) : (23) 
A 2 
] 73C 


The resonant susceptibility may be expressed in terms of its 


peak value, Xp! by: 


ae 


es 
XR Ze = Ay ae) 


where Aw = Wy) 7W57W 


For weak signals in which Xo <<Xyr’ the absolute value 


squared of ye 1s given by 
2 eo 2XpyAw 
e) P P 
Se ee ee (25) 
NR" 2Aw-iy NR Mic: 


oo 


The magnitudes of y at its extrema are XureXp 2° Bete 


Signal strength proportional to the difference in the peak and 
background susceptibility 


2 


Xp 2 
+ >) = Xyp + XypXp 


ANR ~ 2 Be 


T 0. ( 


aus, for Xp ae Xr’ the anti-Stokes power above that due 


to Xyp alone, is 


24 





z 2 
ee ne pho lp ae 


where 1n tiseatmertrterlency cOoerricient to take into account 


a 
the optical losses, detector inefficiencies and observed con- 
version efficiency factor. 

Typically, the largest source of noise when performing 
CARS on samples at one atmosphere pressure is associated with 
the power fluctuations of the pump and Stokes beams. If the 


power of the laser source, P has a relative standard devia- 


ig 
Gaon Of Eo! the standard deviation of the power is Sr The 
fluctuation of the background signal power Pur due to Xyr 
is thus: 
2 2 
aK 3 = K Xup P, Pon, (28) 


assuming each beam has the same relative standard deviation ce 
For a low intensity CARS signal in the absence of back- 
ground one finds the major source of noise to be the shot noise 
of the signal. For an N_ photon signal the standard deviation 
is VN. 
For a pulsed signal the noise power due to shot noise is 


Phy 
(a: 
Py = = (29) 





The signal-to-noise ratio for the background limited case 


given by equations (27) and (28) becomes 


Xp 


on 


S/N 
p*NR 


(30) 


whereas the signal-to-noise ratio for a situation in which 


there is only a very weak signal is given by equations (27) 


Pe) 





and (29) to be 


NS) 
hol 


Sy Nee =k Xp P Z 


l Nm? Ce7 hw). (31) 


As an application of the signal-to-noise problem consider 


the measurement of Xp from a spectral envelope. The uncer- 
Earnty of Xpe Ons 1s expressed as: 2 2 

Oo o 

P 1 
where oF 1s the standard deviation in determining 5 from 


the a pulse of a pulsed laser. For CARS with a background 
limited situation and Xp *<Xyp the uncertainty in the signal 
amplitude is approximately a constant, oer associated with 


the laser power fluctuations. 
3X 
Thus De: = 0 a where ox 1s basically the 


m 
or 
(w, J,T) 


noise power given by equation (28): 


= Z 2 
gee 3 - K Xyr Pa P.-: (33) 
* (ws J,T) 1s the intensity of the spectral envelope at Wes 


due to the quantum rotational number J, at temperature T. 


The result of this analysis yields [16]: 
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In general any standard deviation is dependent on the number 


of times, N, the measurement is taken such that 


i 
on 
PON 
. ° ° LED 6 ° 

The conclusion is that for a given S/N the term (a5) 1s the 
effective number of pulses required at intervals of Aw to pro- 
duce an uncertainty, Ta in the measurement of Xp- 

The uncertainty in determining the rotational temperature 
1s a product of the nearly constant background uncertainty, 


On! and the small change of intensity with temperature. A 


convenient definition for On 1s: 


o. =I 6 (36) 


where 6 is the relative uncertainty in the amplitude of 
max 
the line of maximum intensity. The uncertainty in the tempera~ 


ture then becomes [16]: 





a1 
el ne 
z= 3 7 ). (37) 


Ak m 


The intensity of a line is proportional to the Boltzmann 
population distribution. The change of intensity with tempera- 


ture has been shown [16] to be: 


PE Oneal cL | 


dT gr, o ep? 





(38) 


The resultant standard deviation in temperature is 
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Oo ee ma (39) 


T : 2 
meg Pes 2) 5 
3)I_ KT 


max 


It has been estimated [16] that for an uncertainty in the 
maximum intensity equal to 2% the temperature of CO at 1500°K 
may be determined to within 10°K. 

| This uncertainty in temperature is for an unrealistic 
Situation in that each rotational line is completely re- 
solved. A more realistic situation is the case of 
a continuous envelope as given in equation (14). A computer 
program (see page 59) was written to calculate the uncertainty 


in a temperature measurement, Ons based on equation (37). 


ae : eet 
= 2 - 2 7 oT 
Ab m 
o =TI1I¢e¢ = K 
m  ~p°p Xp “p 
5 2 
i Men. ; = o(K x 07) | 
= 2 : Z 1 I, oT 
ak —p 
¢ 2 
i ee! r X (ws a) 
e 2 2 1]o0 =a oT 
a p Xp 


For the same conditions used to calculate the previous 
uncertainty in temperature the results of the envelope calcu- 
lation yields an uncertainty of 15°K. This value is of the 


Same order of magnitude as the statistical prediction of 10°K, 
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estimated by using equation (39). It is to be expected that 
the uncertainty would be higher for measurements on the 

spectral envelope. Verification of the statistical approach 
means that the much simpler method for predicting the uncer- 


tainty 1s a good approximation. 
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‘V. CARS IN A PLASMA 


LOw pressure ionized gases, or plasmas, represent a class 
of materials that are receiving very close investigation. The 
population of the various energy levels reveals much about the 
plasma. Determining the population of levels is a very diffi- 
cult task with conventional Raman spectroscopy. CARS has been 
utilized to obtain spectra of several plasmas at the Naval Re- 
search Laboratory, Washington, D. C. The experience there 
indicates that the technique is not difficult and produces 
useful information. 

Of importance when analyzing CARS in a plasma is how much 
Saemaiuaii power laser beams may be expected to perturb the 
system under investigation. The mer orescurce of interaction 
between the plasma and laser beams appears to be from the in- 
verse bremsstrahlung process [23]. The free electrons in the 
plasma will gain kinetic energy by this process. 

If the laser frequency is above the plasma frequency the 
light will propagate through but there will be some attenuation 


Or absorption. The plasma frequency is given by [24] 


L 
ame 2 
WO, = oe (40) 
e 
where is 1s the number density of electrons ee e 


1s the charge of an electron (esu) and M, is the mass of an 
electron (grams). 
An estimate of the change in the electron kinetic energy 


may be made by assuming there is an absorption coefficient, a, 


SU 





for the plasma. The energy absorbed per unit volume of the 
plasma is 


a 
Io (l - e ) AtA PjAta 


si) a. A (41) 


where I, is the incident intensity, Py 1s the incident 
power, At is the pulse duration, A is the beam cross sec- 
tion and x is the path length of the medium. If all the 
absorbed energy is transferred to the kinetic energy of the 


electrons then the energy absorbed per electron is 


aP At 
Oo 


EABS 7 An, ' ear 





As a fraction of the initial energy, 3/2kT, 


EABS : aP At ae 
3/2kT 3/2An kT, 


where Te is the electron temperature. The absorption co- 


efficient, a, for inverse bremsstrahlung has been reported 


as [25] 


aL ZG 
D ee . 
3/23 


hv 
fl - exp (- =>) ] (44) 
he M, kT 


re 
tee arr 
where n, is the ion number density, z is the charge of the 
ions and v is the frequency of the incident light. Evalua- 
tion of equation (43) with a defined by equation (44) results 


in a relative change of temperature equal to 


oT 1.3x10°°n,2°P At (lL - exp (~hv/kT) ] 


—_ = ee (45)} 
ne va (kT) 7/7 
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This agrees favorably with an equation for the change of 


electron temperature given by [23] 


ZS 


9 T 
AkT ax 0 nz [l - exp (“hv/kT) lf Pidt 


—— (46) 
2 


va) 9/4 A 
where (kT) is in electron volts and PS is the power of 
the laser in watts. 

The equipartion time, tee is the average time required 
for the electrons to transfer their excess kinetic energy to 
the ions and thus reach an equilibrium temperature. The equa- 


eran £LOr te 1s given by [25] 


25.2 a(t _)?/2 
ee ne (47) 


Z 
Ne 


where A is the atomic weight of the ions, T. is the electron 


temperature in degrees Kelvin, and Ne is in em >. If the 


duration of the laser pulse is shorter than te then any 
energy the electrons gained by inverse bremsstrahlung will 
not be transferred to the ions and there will be very little 
perturbation of the molecules under investigation. 


Consider the following as a typical example: 


plasma: Deuterium at 50 Torr,Te = 10° K 


7a = Moe aa, Z0 em slong 


laser: 5500A, 5MW, Gmeacerna beam area, 20 nsec 


aduration 


The plasma frequency, Wr equals 1.8 x ice sec and eae 


The equipartion time, toe equals 20 micro seconds. The 


32 





relative change in electron temperature given by equation (46) 
is .0001%. Little perturbation of the plasma is expected from 


these conditions. 
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VI. BACKGROUND SUPPRESSION EXPERIMENT 


As discussed earlier, the nonresonant susceptibility, Xp’ 
is responsible for the background signal. The power amplitude 
fluctuations of the lasers produce noise on this signal. The 
resonant Signal above the background also has noise on it. It 
would be much easier to detect the resonant signal if there 
were no noise on the background or, better yet, no background 
at all. The desire for reduction of the background signal led 
to an idea [26] for an experiment to determine if it was possi- 
ble to completely suppress the background. In addition to 
yielding much better signal quality the elimination of the 
background would limit the detectibility of CARS to the limit 
of conversion efficiency. 

Typically, the diluent of a sample provides the largest 


Semeribution to The basic idea for suppression of the 


XNR’ 


background signal was to perform two CARS experiments in series. 
In the first experiment the generated anti-Stokes beam, Was 
would be due to only XnR (On a calvcnt mand im tne Cteher 


t 
experiment the anti-Stokes beam, War would be due to cf 


XNR 


(the sample plus the diluent). If the beam, from the 


a3 


First cell, were then forced out of phase by Na radians 


XR 


t 
(where N is an odd integer) relative to the second beam, Wor 
and superimposed on the second then the resultant signal should 


be only that due to of the sample. 


XR 
Molectron Corporation of Sunnyvale, California, through 


Dr. G. Klauminzer provided equipment, facilities and assistance 


for conducting the experiments. 
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Figure 2. Arrangement of Ban temeae for Background Suppression 
Experiment. 


Figure 2 shows the basic arrangement found to be most 
convenient. A pulsed nitrogen laser was used to drive two 
tunable dye lasers in an arrangement reported by Dr. 
Klauminzer in Reference 12. 

Benzene was chosen as the first material to work with 
Since it was known that the CARS signal was very strong and 


W. and w were 


visible. For benzene the values for w 9 3 


‘ig 
480., 504. and 458.2 nm, respectively. 


Both cells were filled with benzene in anticipation of 


demonstrating complete cancellation of W with We 


known that a translation of cell one along the focal axis 


It was 


would change the intensity of W by changing the volume in 
the focal region but it was also discovered that the same 
LO 


' 
translation changed the phase of relative to 


Ww Wa 
f 
completely cancel W it was known that several conditions 
had to be met. The beam had to be tm radians out of phase, the 


intensities had to be equal and the beams had to be coinci- 


dent. A plot of the resultant intensity versus the translation 


25 


oe 








of cell one revealed that there was constructive and destruc- 
tive interference of the two beams. The destructive interfer- 


ence was not complete though. A visual observation of the 


(On ot) 
3 


g 
W. beam off of the We 


changed. It was suspected that chromatic aberration in the 


beam revealed interference fringes. By steering the 


beam the fringe pattern could be 


focal lens was responsible. Calculation of the effect of 
chromatic aberration indicated that the fringe pattern should 
be approximately that observed. The use of an achromat for 
the focal lens reduced the number of fringes visible but it 
was found that a slight rotation of the lens changed the 
pattern considerably. 

The best measured ratio of minimum to maximum intensity 
as a function of cell translation was observed to be 10:1. 
Since it was the electric fields that were interfering, the 
10:1 ratio amounts to about a 60% cancellation factor. The 
cancellation factor was arrived at by considering the incom- 
plete phase cancellation of the electric fields. The observed 


intensity 1S proportional to 
E, e + E, e . (48) 


If the amplitudes Ey and E. are equal, the resulting 


intensity 1S proportional to 


ie 46. ene ae 
Eeierre ule. ¢ +e 2) (49) 


OF 


Ia E*(2 + 2 cos(0, - 0,)). (50) 


2 
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Complete constructive interference, when 0, = Q results 


al ye 


in an intensity four times that of one beam, and complete de- 
structive interference would result in no intensity. 
It is not surprising that on the average complete inter- 


ference does not take place. For a small difference in 04 


and 0. the observed magnitude is slightly less than 4E7 or 


slightly more than zero. 


An observed 9:1 intensity ratio is equivalent to a 3.6 to 


.4 ratio on the normalized scale, ne = 1. Thus, the inter- 


ference produces a signal equal to 40% of what one beam would 
produce, therefore, a cancellation factor of 60%. 

Next, a three molar sodium benzoate solution was used in 
the cells in place of benzene. Again there was incomplete 
cancellation. The best observed ratio of intensities was 5:l. 


It had been predicted that sweeping the frequency of Wo to 


record a spectrum would not alter the phase difference of W 2 


| e 
and Wa The first observations made while sweeping w re- 


2 


went in and out of phase over a 


vealed that Ww and W 
period of about 10 cmt. By changing the optical path length 
Ot: Ls Og Wy and Wo by use of microscope slides the 
period could be increased to about 60 cm > 

The observed increase in the period of the oscillation 
verified that the optical path length could be controlled to 
shift Ww. and Ww, in and out of phase. With the appropri- 
ate optical path length the period of the oscillation could 
be adjusted to infinity. 


A spectrum was recorded for the three molar sodium benzoate, 


two cell, series cancellation experiment. Of interest from the 


oy 





spectra is the interference of the background with respect to 
the noise levels observed. Figure 3 is a portion of that 
spectrum depicting the observed resultant signal. When the 
interference was constructive the noise was much higher than 
for the destructive interference. 

As stated earlier, the reason for the experiment was to 
increase the sensitivity of CARS. An analysis of the signal- 
to-noise ratio for a cancellation experiment indicates that 
the cancellation has to be much more effective than was ob- 
served in this experiment to gain an order of magnitude in- 
crease in the signal-to-noise ratio. 

Assume oe is the fraction of background intensity remain- 


ing after cancellation. Then, the S/N ratio is given by 


/N " 
S/N * sey 
3€, (Xu 


Similar to the development of equation (30). An increase in 
S/N of 10 requires a = .01l, or a 99% cancellation factor. 
A 400:1 intensity ratio would correspond to such a cancella- 
EPoneractor. 

The experiment as performed was very sensitive to the 
alignment and quality of the optics involved. One of the 
major contributions to the incomplete cancellation is probably 


the diffraction limited ability of the center lens to perfect- 


iar 


factors such as turbulence and inhomogeneity of the samples 


ly image the initial W beam onto Other macroscopic 


are other obstacles to complete cancellation. 
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Vet eereoouls AND DISCUSSTON 


The basis for determining the areas of application for a 
new analytical tool, such as CARS, must depend on recognizing 
the limits beyond which the technique cannot be expected to 
yield the desired information. 

The results of this analysis are presented as lower limits 
of concentration that are expected to be detectable by CARS. 
For a mixture of gases, the detectability of one gas depends 
on the ratio of Xp and Xn! the same as Taran's plot for 
H. and N.- Plots of the ratios for other gases with N, as 
a diluent are presented in Figure 4. The minimum detectable 
concentration is dependent on the required signal-to-noise 
ratio for a given experimental condition. The lower limit of 
detection for a pure gas at reduced pressures would be limited 
to the conversion efficiency if not for the photon statistics 
involved. From the signal-to-noise analysis, estimates of the 
feeectavle lower pressures are summarized in Figures 5 and 6. 

Based on the results of Section IV the rotational tempera- 
ture of a gas is expected to be measurable to within 1%. That 
assumes the amplitude fluctuations of the CARS signal observed 
remain less than 10%. 

In a plasma, the inverse bremsstrahlung process 1s ex- 
pected to be a relatively insignificant perturbation to the 
parameters determined by a CARS experiment. 

The experimental attempt at cancellation of the background 


confirmed the idea that interferometric cancellation is a 
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possible technique for improving CARS when the background be- 
comes a problem. The observed results were not of the magni- 
tude hoped for but did serve to verify the concept and to 
recognize a possible source of limiting performance. 
Turbulence in the samples and inhomogeneity may be part 

of the reason that incomplete cancellation occurs. The major 
contribution to the lack of complete cancellation is probably 
the result of diffraction limited optics. The first beam 
cannot be refocused through a lens and have precisely the same 


geometry that it had originally. 


4} 





Relative Anti-Stokes Intensity 


Figure 4. 


in’ 


Concentration [ppm] of Sample 


Log-log plots of ge gas mixtures at 1 atmosphere 
total pressure. at 300 K. Upper CO plot at 


300 K, Lower me: Ai me OO i Ke Xue or N. assumed 


Pepeieoee on. C8] at 300K, and 
“-18 ene mee 
eo a xO erg at 1500 XK. 
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SIGNAL-TO-NOISE RATIO 


Peragure 5. 


I es 
LYS 


| 
KA 
in? n° 104 io 102 
SAMPLE PRESSURE (ATMOSPHERES) 


Signal-to-noise ratios for gas mixtures at 300 K 
with a total pressure of 1 atmosphere. Background 
pressure and signal due to N.. The area boundar- 
ies are defined by assuming values of the un- 


certainty in the pulsed laser intensity, E from 


eo on. 
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Signal-to-Noise Ratio 


Figure 6. 


Ho 0» C0, No 





a i n! 


10 10 


PRESSURE (ATMOSPHERES) 


Signal-to-noise ratios for Doppler broadened 

gas at 300°K as a function of pressure. Noise 
due to photon statistics. The assumed parameters 
are 
Np 
X = 5500 A. 


= .0l, Py = 5MW, P. = .5MW, T = 20 nsec, 


a4 








Verte CONGCEUS IONS 


This analysis presents a model by which the limits of the 
information available from CARS may be predicted. At STP 
typical diatomic gases should have a signal-to-noise ratio of 
l at partial pressures of about om atmosphere. Hydrogen, be- 
cause of its narrow linewidth and population distribution, 
should have a signal-to-noise ratio of l at about ome atmos- 
phere. Pure gases at low pressure where the linewidth is 
determined primarily by Doppler broadening, may be detected 
at pressures several orders of magnitude lower than this. 

The ultimate detectability depends upon laser power and satura- 
tion effects. 

Simple relationships have been derived which allow the 
prediction of the uncertainty in temperature of a gas for an 
assumed signal-to-noise ratio. 

In a plasma the inverse bremsstrahlung process is found 
to transfer insufficient energy to the plasma to significant- 
ly modify the parameters measured. 

The experimental attempt at background reduction by inter- 
ferometric cancellation verified that the principles were 
correct. Complete cancellation, however, was not observed. 
This is thought to be due to the diffraction limitations of 
optical lenses. 

One area of additional research and experimentation that 
could further extend the applications of CARS is that of the 


background suppression. Another area of improvement needs to 
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be that of the fluctuating intensities of the high power 


pulsed lasers that are available. 
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APPENDIX A 
STIMULATED RAMAN SCATTERING 


A. DEVELOPMENT 

The interaction of two electromagnetic fields in a material 
is such that different combinations of the two may produce 
additional fields. The two additional fields of interest here 
are the Stokes and anti-Stokes frequencies. 

Raman scattering has been observed extensively and found 
to produce the Stokes and anti-Stokes frequencies at plus and 


Minus the vibrational frequency, w of the molecule from 


V7 
the incident light, w 


I 
stokes “1” “v (Al) 
Manti-Stokes “I = aay (A2) 


Stimulated Raman scattering predicts a gain in the strength 
of a field as it traverses a media. The media of interest in 
this study are gases. The gain will be related to the non- 
linear susceptibility of those gases. 

Beginning with Maxwell's equations, which contain various 
assumptions implied by their form, the plane wave equation is 
derived from the curl of the curl of the electric field, E. 
Assuming a solution to the wave equation of the form: 

l -i(kz - wt) 
E=5 (ee aC Ge) (A3) 
where e« is the amplitude, k is the propagation constant 
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and w is the frequency. For simplicity, all fields have 
been reduced to one dimension along the 2 axis. From the 
nonlinear polarization term the various possible frequencies 
are defined and their field properties (gain or attenuation) 
are explained by computing the solution to the wave equation. 
For completeness both MKSA and CGS results are presented. 

The development in MKSA units begins with Maxwell's basic 


and constituent equations. 


Vos o= Jot i ee 
VxE= -B J = oF 

Ve E= 0 EY = eax Bor Put 
Ve B= 0 B= uh 


Vx H=ocE+ ec E+tP = cE + cE + Pi. 
O NL 


Then the wave equation, defined from VX(VxE) is 


2 d eee 
oe = 0 A US ae (VY x H) 
Or 
2 “s 
= gt ee 
Es | eer re 
O 
Using the assumed form of the solution, equation (A3), and 


2 , 
the assumptions that o= 0 and 39 6/927 is very small, 


then equation (A4) becomes 
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mye co ee aC, tf . = (cau ezee oe ene.) 
2n 
ke“ emer 2 = wt) Dee 
- ae 1 De < qf eG Ae == a PL’ (A5) 
The usual definitions: 
2 1 
c= : 
Koo 
n = ABS 
k = w n/fe 
have also been used. 
For the two interactive fields Ey and E, the nonlinear 
porarization, Por, is expressed as: [18] 
a (3) 2 
Three frequencies are possible from the EE, combination; 
Was 2W1—W572u,+uQ,- It ES is the Stokes field with frequency 
Wo (assuming w4< We) then a resultant field, Ea, of fre- 
quency 204 -Wo is the anti-Stokes field. For the anti-Stokes 
calculations Pur is given by: 
-i[((2k,-k.,) 2— (2w,-w..) t] 
e By al: * a2 1 ~2 
PML = 3E OX (SEE LE> [e Gare Jol. (A7) 


The resultant field may be ES for the case when 


“i (k,2 - Wot) 


_ Bik x 
Pur = 3€ UX (SE, £4 e,) le GG cli (A8) 
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1. Stokes Gain Equation 


To solve equation (A5) for the Stokes condition, equa- 
tion (A8), there are two possible alternatives that lead to a 
gain equation. One way is to assume that 9€,/ 92 is very 


small and Ky is complex such that 


— I : 7? 
Kk, Sec 1K 


The imaginary part of k k", is then evaluated as: 


Dt 


wy 3x"|E,°|? 
a) Z 1 
k = 
Zen 


where yx" is the imaginary part of (3) and E, is the 


magnitude of Ey at z= 0. 


The observed intensity of E is proportional to the 


2 
time average of jE, |? or 


== 1._* 
JE,|" = 5EE 
where 
a * i ik’'2 -k"z Ror eZ eZ 
5E5E, = g lene ce + E, e 2 aCe.) 
=— -2k"Z re 
aS ; JE. | is proportional to e . x" is a negative 
Guantity [27] and therefore ance 2 is an exponentially 


Gaovengetunetron with distance. 
The other possible way to solve equation (A5) is to 


let 


and then solve for J€5/ 92. The result is that 
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J€ W 2 
2” 2 : ae. Oo 2 
fee EX) IED I Fe! 


and again for yx" negative there is gain for increasing 2. 
Both approaches result in an intensity gain coeffici- 


ent, g, such that 





me Oo. gz 
I. I, = 
where Ls is the intensity of I, ac. Za 05, rhe. result 1s 
phat 
05 S Z 
5. 3 Toe [3x"| |B, | . (A9) 
For cgs units the gain is 
4tw 2 
2 , O 
Seeteeeo ok | Ey |. (A10) 





es 
2. Anti>Sstokes Gain Equation 
Solving equation (A5) for the anti-Stokes condition 


(A7) results in the gain equation for the anti-Stokes field. 


The only way to solve equation (A5) is to assume 


. 25 2 
eS Gas 
nna ae 
c 
The result is that 
J€ 1w 
Se ee * iAkz 
ee “ W5e 3(x' + ix") 38121 @& (All) 
where Ak = 2k, = Ky - K3- 
Again if y" is negative, there is gain. In cgs units 


Dek 





J€ J€ 


3 ye 3 
dz(cgs) — aos OZ ) KSA. (Al2) 


B. CORRELATION OF VARIOUS ARTICLES 

The only published values for nonresonant susceptibilities 
of gases are those given by Rado [21]. It was desirable to 
check for consistency in the equations used by Rado and the 
ones used to calculate resonant susceptibilities. 


The basis for Rado's experimental values of Xyp Was the 


resonant value of re OIG H. from other experimental 
data [23]. The gain coefficient for Stokes Stimulated Raman 
Scattering provided Rado with the needed value of Xp fox H. 


at pressures greater than 10 atmospheres. The gain coeffici- 


ent was given in [28] in cgs units as 


Kk (8n°L A No“ } 5 


wmere iL is the local field correction factor, N is the 
number density of atoms, A is the population ratio, i is 
the phonon lifetime, a is the polarization coefficient and 
Ty is the incident laser intensity. 

An independent article by Maier [29] used the following 


definition for the gain coefficient 


scien SmN ne 
ot rie = Av 


where Av is a linewidth. Rado's gain is equivalent to Maier's 


lek 





and 


a 


e iene Av. 


The experimentally observed gain for Ho above 10 atmos- 
pheres was stated as 1.45 x ome em/W for an incident light 
source of 975 nm. If the gain coefficient developed in equa- 


tion (Al10) is divided through by intensity then 





ZZ 
Pie uses | 
BOS, wee 2 
= 75 ; (A15) 
iL Crm 
Solving for |3x"| the result is 
ger 
i" 2 
ae ee (A16) 
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assuming ne = 1. The value for 3x" is the same as that 
used by Rado to represent yx". 

The gain coefficients given in two textbooks more nearly 
resemble the gain coefficient in equation (A9). Pantell and 


Putoff [30] define g in MKS units as 


4n x" iT 
SS err (Aa) 
pu has 


whereas equation (9) may be rearranged as 


4m |3x"| T, 


g = ———_—— (A18) 
AoNoEQS 


differing from Pantell by a factor of Ny: Hier ny, = 1, they 
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become the same. Yariv [27] gives a voltage gain coefficient 





as 
k 
' Z 8 
gy lel” Ix"). (A19) 
ND, 
The intensity gain, g, would be 2g'. Thus, the intensity 
gain coefficient is 
k WwW 
2 2 #9 2 Ty Z 
oe Se a EL Ke le ee a NE (A20) 
N5 2 
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FROGRAM TO CALCULATE THE LOCAL FEAK VALUES OF TFE 
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